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Abstract 

Detailed mean field and Monte Carlo studies of the dynamic magnetization- 
reversal transition in the Ising model in its ordered phase under a com- 
peting external magnetic field of finite duration have been presented here. 
Approximate analytical treatment of the mean field equations of motion 
shows the existence of diverging length and time scales across this dy- 
namic transition phase boundary. These are also supported by numerical 
solutions of the complete mean field equations of motion and the Monte 
Carlo study of the system evolving under Glauber dynamics in both two 
and three dimensions. Classical nucleation theory predicts different mech- 
anisms of domain growth in two regimes marked by the strength of the 
external field, and the nature of the Monte Carlo phase boundary can be 
comprehended satisfactorily using the theory. The order of the transition 
changes from a continuous to a discontinuous one as one crosses over from 
coalescence regime (stronger field) to nucleation regime (weaker field). Fi- 
nite size scaling theory can be applied in the coalescence regime, where the 
best fit estimates of the critical exponents are obtained for two and three 
dimensions. 
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1 Introduction 



The study of the response of pure Ising systems under the action of a time- 
dependent external magnetic field has been of recent interest in statistical physics 
[000. A whole class of dynamic phase transitions emerged from the study of 
such driven spin systems under different time dependences of the driving field. A 
mean field study was initially proposed by Tome and Oliveira [|J where the time 
dependence of the external perturbation was periodic. Subsequently, through 
extensive Monte Carlo studies, the existence of a dynamic phase transition under 
periodic magnetic field was established and properly characterized @[f| 0. Later, 
efforts were made to investigate the response of such systems under magnetic 
fields which are of the form of a 'pulse' or in other words applied for a finite 
duration of time. All the studies with pulsed fields were made on a system 
below its static critical temperature T c °, where the equilibrium state has got a 
prevalent order along a particular direction. The pulse is called 'positive' when it 
is applied along the direction of the prevalent order and 'negative' when applied 
in opposition. The results for the positive pulse case was analyzed by extending 
appropriately the finite size scaling technique to this finite time window case, and 
it did not involve any new phase transition or introduced any new thermodynamic 
scale || . However a negative field competes with the existing order and depending 
on the strength h p and duration At of the pulse, the system may show a transition 
from one ordered state with equilibrium magnetization +mo (say) to the other 
equivalent ordered state with equilibrium magnetization — m ||. This transition 
is called here the "magnetization-reversal" transition. It may be noted that a 
magnetization-reversal phenomenon trivially occurs in the limit At — > oo for 
any non vanishing value of h p at any T < T 6 °. However, this is a limiting case 
of the transition, which is studied here only for finite At. In our studies the 
magnetization-reversal need not occur during the presence of the external field. 
In fact, it will be shown later that the closer one approaches to the threshold 
value of the pulse strength longer is the time taken by the system, after the 
field is withdrawn, to relax to the final ordered state. We report here in details 
the various results obtained for this dynamic magnetization-reversal transition in 
pure Ising model in two and three dimensions. 

The model we studied here is the Ising model with nearest neighbour inter- 
action under a time dependent external magnetic field, described by the Hamil- 
tonian 

H = -\ E J ii 3 A - E Ht)Si, (i) 

where is the cooperative interaction between the spins at sites % and j respec- 
tively and each nearest-neighbour pair denoted by [. . .] is counted twice in the 
summation. We consider the system at temperatures only below its static critical 
temperature (T < T c °). The external field is applied after the system is brought 
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to equilibrium characterized by an equilibrium magnetization ttiq(T). The field 
is uniform in space (hi(t) = h(t) for all i) and its time dependence is given by 



h(t) = -h p , for t < t < t + At (2) 
= , otherwise. 

Typical responses of the time dependent magnetization m(t) under different h{t) 
are shown in figure 1. As mentioned before, for appropriate combination of 
h p and At, magnetization-reversal transition occurs when the system makes a 
transition from one ordered state to another. This transition can be observed 
at any dimension d greater than unity for systems with short range interactions. 
This is because one has to work at temperatures T < T c ° where, in absence of a 
symmetry breaking field, the free energy landscape has got two equivalent minima 
at magnetizations m = ±m . A phase boundary in the h p — At plane gives the 
minimal combination of the two parameters at a particular temperature T(< T°) 
required to bring about the transition. 

A full numerical solution as well as an analytical treatment in the linear limit 
of the dynamic mean field equation of motion shows the existence of length and 
time scale divergences at the transition phase boundary The divergence of 
length and time scales is also observed in Monte Carlo (MC) simulation study of 
Ising model with nearest neighbour interaction evolving under a negative pulse 
through single spin flip Glauber dynamics The phase diagram for the tran- 
sition was obtained for both MF and MC studies. While the phase boundaries 
for the two cases are qualitatively of similar nature, there exists a major differ- 
ence which can be accounted for by considering the presence of fluctuations in 
the simulations. In the MC study, there exists two distinct time scales in the 
problem : (i) the nucleation time tn is the time taken by the system to leave 
the metastable state under the influence of the external magnetic field and (ii) 
the relaxation time tr is the time taken by the system to reach the final equi- 
librium state after the external field is withdrawn. While is controlled by 
the strength h p of the external pulse and is bounded by its duration At which 
is finite, tr is the time scale that diverges at the magnetization-reversal phase 
boundary. According to the classical nucleation theory (CNT) [[□]], there can be 
two distinct mechanisms for the growth of domains or droplets depending on the 
strength of the external field. Under the influence of weaker external magnetic 
fields, only a single droplet grows to span the entire system and this is called the 
single-droplet (SD) or the nucleation regime. On the other hand, under stronger 
magnetic fields, many small droplets can grow simultaneously and eventually co- 
alesce to form a system spanning droplet. This is called the multi-droplet (MD) 
or the coalescence regime. The crossover from SD to MD regime takes place at 
the dynamic spinodal field or husp{L,T) which is a function of system size L 
and temperature T. The nucleation time changes abruptly as one crosses over 
from SD to MD regime even along the same phase boundary. The nature of the 
transition too changes from a continuous one in the MD regime to a discontinuous 
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nature in the SD regime. All our simulation observations for the dynamic phase 
boundary compare well with those suggested by the CNT. The investigations 
about the relaxation time tr and the correlation length £ are also discussed here. 
The application of scaling theory in the MD regime gives the estimates of the 
critical exponents for this dynamic transition. The organization of the paper in 
as follows : We discuss the MF results in the next section and the MC results for 
square and simple cubic lattices in section 3. A brief summary and concluding 
remarks are given in section 4. 

2 Mean field study 

The master equation for a system of N Ising spins in contact with a heat bath 



evolving under Glauber single spin flip dynamics can be written as \T2 

jP(S u ---,S N ;t) = -J2W j (S j )P(S l ,...,S N ;t) 

+ J2W j (-S j )P(S 1 ,... 1 -S j ,...,S N ;t), (3) 
j 

where P (Si, ■ ■ ■ , Sn] t) is the probability to find the spins in the configuration 
(Si, ■ ■ , Sn) at time t and Wj(Sj) is the probability of flipping of the jth spin. 
Satisfying the condition of detailed balance one can write the transition proba- 
bility as 



was,) = i 



1 — Sj tanh 



'Ei JijSi(t) + h 



j 



T 



(4) 



where A is a temperature dependent constant. Defining the spin expectation 
value as 

m i = (Si) = ^2S i P(S 1 ,...,S N ;t), (5) 

{S} 

where the summation is carried over all possible spin configurations, one can 
write 

A^i = - m , . + (tanh j&M^j) . (6) 
Under the mean field approximation (|J) can be written after a Fourier transform 



as 



X^- = -m q (t) + tanh M*Wg + M*A ? (7) 
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where J(q) is the Fourier transform of J^. Equation (|7|) is not analytically 
tractable and one can only look for solutions in the small m q limit where terms 
linear in m q are dominant. The linearized equation of motion, therefore, can be 
written as 



dm q (t) 
dt 



(K(q)-l)m q (t) 



T 



(8) 



where K(q) = J(q)/T. When we are concerned only with the homogeneous 
magnetization, we consider the q = mode of the equation and writing m q=0 = m 



and h 



q=0 



h, we get 



dm 
~dt 



(K(0)-l)m(t) + 



T 



(9) 



In the mean field approximation K(0) = T C MF /T with T C MF = J(0) and for small 
q, K(q) ~ -^(0) (1 — q 2 ). Differentiating (0) with respect to the external field, we 
get the rate equation for the dynamic susceptibility Xqif) as 



A 



dx q {t) 



-x q (t) + 



f J(q)x q (t) + 1 



dt vyv \ T J 

which in the linear limit can be written as 



sech 



J(q)m q (t) + h q (t) 
T 



(10) 



dt 



X 



(K(q) - 1) Xq (t) + 



T 



(11) 



Before we proceed with the solutions of these dynamical equations, we divide 
the entire time zone in three different regimes : (I) < t < to, where h(t) = (II) 
to <t < to + At, where h(t) = -h p and (III) t + At < t < oo, where h(t) = 
again. We note that ([5]) can be readily solved separately for the three regions 
as the boundary conditions are exactly known. In region I, dm/dt = and the 
solution of the linearized (|9p becomes trivial. We, therefore, use the solution of 
(0) in region I (mo = tanh {rriQT^ IF /Tj) as the initial value of m for region II. 
Integrating M) in region II, we then get 



m(t) 



hp 
AT 



+ m 



h p s 
AT 



exp [bAT (t - t )] 



(12) 



where b = 1/AT and AT = T C MF — T. It is to be noted that in order to justify the 
validity of the linearization of (^) one must keep the factor inside the exponential 
of ( |12"D small. This restricts the linear theory to be valid at temperatures close 
to T C MF and for small values of At. Writing m w = m(t + At), we get from ( |12"D 



hp 
AT 



+ [mo- 



il 



AT 



!■ 1 e bATAt 



(13) 
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It is to be noted here that in absence of fluctuations, the sign of m w (h p , At) solely 
decides which of the two final equilibrium states will be chosen by the system 
after the withdrawal of the pulse. At t — t + At, if m w > 0, the system goes back 
to +m state and if m w < 0, magnetization-reversal transition occurs and the 
system eventually chooses the — m state (see figure 1). Thus setting m w = 0, we 
obtain the threshold value of the pulse strength at the mean field phase boundary 
for this dynamic phase transition. At any T, combinations of h p and At below the 
phase boundary cannot induce the magnetization-reversal transition, while those 
above it can induce the transition. From ( |13|) therefore we can write the equation 
of the mean field phase boundary for the magnetization-reversal transition as 



K(m,t)= i ^_Zm - (") 



Figure 2 shows phase boundaries at different T obtained from (|T4l) and compares 
those to the phase boundaries obtained from the numerical solution of the full 
dynamical equation (^) . The phase boundaries obtained under linear approxima- 
tion match quite well with those obtained numerically for small values of At and 
at temperatures close to T C MF , which is the domain of validity of the linearized 
theory as discussed before. In region III, we again have h(t) = and solution of 
(H) leads to 

m{t) = m w exp [bAT {t - (t + At)}]. (15) 

We define the relaxation time Tr F , measured from t = to + At, as the time 
required to reach the final equilibrium state characterized by magnetization ±mo 
in region III (see figure 1). From (|I3|) therefore we can write 



MF 




J^MF _ rp 



In \m w \ . (16) 



A point to note is that m(t) in ([15]) grows exponentially with t and therefore 
in order to confine ourselves to the linear regime of m(t), mo must be small (T 

close to T C MF ) and t < r^ IF . The factor (t^ if — T^j gives the usual critical 
slowing down for the static transition at T = T^ IF . However, even for T -C T^ F , 
Tft IF diverges at the magnetization-reversal phase boundary where m w vanishes. 
Figure 3 shows the divergence of t^ f against m w as obtained from the numerical 
solution of the full mean field equation of motion (^) and compares it with that 
obtained from (|16|). 

Solution of Xq{t) is more difficult as all the boundary conditions are not di- 
rectly known. However, Xg{t) can be expressed in terms of m(t) and the solution 
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of the resulting equation will then have the t dependence coming through m(t), 
which we have solved already. Dividing ([T^) by (0) we get 



(17) 



dm(t) 

which can be rewritten in the linear limit as 



Xq (t) + ( J ^f +1 ) 


sech 2 


' J(q)mq(t)+h q {t)~ 
T 


m 9 (t) +tanh( J(,?)m « 


t)+h q (t)\ 





dXq 



dm 



l/T(K(q) 



1) and a„ 



q m + h{t)/AV 
(K(q) - 1) / (K(0) 



where V 
small q. 

In region II, solution of ( |l8l) can be written as 



Xq (t) = -r + ( Xq + r) 



m(t) - h p / AT 
itiq — hp/ AT 



q 2 / AT for 



(19) 



where x\ is the equilibrium value of susceptibility in region I. Solving ([Tj|) in 
region III with the initial boundary condition m (to + At) = m w , we get 



x q {t) 



-r + ( x , (to + At) + r) 

m(t) V 



r+ xS + r 



/m(t)\^ 



bATAta q 



(20) 



where use has been made of ([19]) and (p~3T) . The dominating g dependence in Xqif) 
is coming from (l/m w ) aq when m w — > as one approaches the phase boundary. 
The singular part of the dynamic susceptibility can then be written as 



X<(t) = Uo + r exp 



_ q2 feMF 



(21) 



where for small values of m w the correlation length £ MF is given by |L0 

1 



i MF = e iF K 



Tc 
AT 



In 



\m,. 



(22) 



Thus the length scale also diverges at the magnetization-reversal phase boundary 
and this can be demonstrated even using the linearized mean field equation of 
motion. Equations (pf ) and (|22|) can now be used to establish the following 
relation between the diverging time and length scales : 



MF T 

Tr ~T, 



(23) 
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which leads to a dynamical critical exponent z = 2. It may be noted that these 
divergences in t^ if and £ MF are shown to occur for any T < T^ IF , and these 
dynamic relaxation time and correlation length defined for the magnetization- 
reversal transition exist only for T < T C MF . 

It may further be noted from ( pip that ^ as £ MF — > oo, thereby 

producing a minimum of x 9 at the phase boundary. The absence of any divergence 
in the susceptibility is due to the fact that at t = t + At, there remains no 
contribution of m w in Xq{t) as is evident from (p0|). However, numerical solution 
of (|i~7|) for g = mode shows a clear singularity in the homogeneous susceptibility 
Xo at the magnetization-reversal phase boundary {m w = 0), as depicted in figure 
4. One can also have a numerical estimate of £ MF by solving (|T7|) for different 
values of g. Figure 5 shows plots of Xqif) against m w for different values of q. 

The inset of figure 5 shows the variation of (^ MF ^ against (In |m„,|) _1 , where 
£ MF was obtained by fitting the data of figure 5 with straight lines. It is clearly 
seen from the inset that for small values of m w the linear approximation agrees 
quite well with the numerical results. 



3 Monte Carlo Study 

We now study the transition using Monte Carlo simulation with single spin- 
flip Glauber dynamics |L3| . Working at a temperature below the static critical 
temperature (T c ° ~ 2.27 and 4.51 [|n]] in units of the nearest neighbour inter- 
action strength J for square and simple cubic lattices respectively), the system 
is prepared by evolving the initial state (say with all spins up) under Glauber 
dynamics for the temperature T. The evolution time to is usually taken to be 
sufficiently larger than the static relaxation time at T to ensure that the system 
reaches an equilibrium state with magnetization mo before the external magnetic 
field is applied at time t = t . The magnetization m(t) starts decreasing from 
its initial value m due to the effect of the competing field during the period 
< t < to + At, and it assumes the value m w at t = t + At. Due to presence 
of fluctuations, m w < does not necessarily lead to a magnetization-reversal 
whereas even for m w > fluctuations can give rise to a magnetization-reversal. 
This is in contrast with the mean field case, where due to the absence of any 
fluctuation the sign of m w solely determines the final state. In the MC study, 
however, on an average the final state is determined by the sign of m w (see fig- 
ure 1). The magnetization-reversal transition phase boundary therefore again 
corresponds to m w = 0. 

Figure 6 shows phase boundaries at different temperatures for square and 
simple cubic lattices. The data points for d = 2 are averaged over 500 different 
Monte Carlo runs (MCR) and those for d = 3 are averaged over 150 MCR. A 
qualitative difference between the MF and the MC phase boundaries may be 
noted here. In the former, even for At — > oo, due to the absence of fluctuations, 
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hp must be greater than the non-zero coercive field to bring about the transition 
and therefore the phase boundaries flatten for larger values of At. However, 
in real systems fluctuations are present and even an infinitesimal strength of the 
pulse, if applied for very long time, can bring about the transition. This is evident 
from the asymptotic nature of the phase boundaries for large values of At. 

It is instructive to look at the classical theory of nucleation to understand 
the nature of the MC phase diagram of the magnetization-reversal transition. A 
typical configuration of a ferromagnet, below its static critical temperature T c °, 
consists of droplets or domains of spins oriented in the same direction, in a sea of 
oppositely oriented spins. According to CNT, the equilibrium number of droplets 
consisting of s spins is given by n s = N exp (— e s /T), where e s is the free energy of 
formation of a droplet containing s spins and N is a normalization constant. In 
presence of a negative external magnetic field h, the free energy can be written as 
e s = —2hs + as^ d ~ l '' d 1 where the shape of the droplet is assumed to be spherical 
and cr(T) is the temperature dependent surface tension. Droplets of size greater 
than a critical value s c are favoured to grow, where s c = [a(d — 1) / (2d \h\)] d is 
obtained by maximizing e s . The number of supercritical droplets is therefore 
given by n Sc = N exp \—A d a d \h\ /T\, where A d is a constant depending on 
dimension only. In the SD regime, where a single supercritical droplet grows 
to engulf the whole system, the nucleation time is inversely proportional to the 
nucleation rate I. According to the Becker-Doring theory, I is proportional to 
n s and therefore one can write 



(x r 



cx exp 



d-l 



[T\h\ 



However, in the MD regime the nucleation mechanism is different and in this 
regime many supercritical droplets grow simultaneously and eventually coalesce 
to create a system spanning droplet. The radius s]J d of a supercritical droplet 
grows linearly with time t and thus s c oc t d . For a steady rate of nucleation, the 
rate of change of magnetization is It d . For a finite change Am of the magnetiza- 
tion during the nucleation time rff D , one can write 



rT MD 

Am oc J N It d dt = I 



d+l 



Therefore, in the MD regime one can write |~5| fl6 



T™ D OC J-W+l) exp 



T(d + 1) \h\ 



During the time to < t < t + At, when the external field remains 'on', the 
only relevant time scale in the system is the nucleation time. The magnetization 
reversal phase boundary gives the threshold value h c p of the pulse strength which, 
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within time At, brings the system from an equilibrium state with magnetization 
+mo to a non-equilibrium state with magnetization m w = 0_, so that eventually 
the system evolves to the equilibrium state with magnetization —m®. The field 
driven nucleation mechanism takes place for to < t < to -\- At and therefore 
equating the above nucleation times with At, one gets the for the magnetization- 
reversal phase boundary 



In (At) = a + C 



h c 
p 

d 



l-d 

, in the SD regime 

(24) 



= c 2 + C hp /(d + 1), in the MD regime 
where C = Ad<j d /T and c\, c 2 are constants. Therefore a plot of ln(At) against 

id-1 

would show two different slopes corresponding to the two regimes [[L7[ . 
? igure 7 shows these plots and it indeed have two distinct slopes for both d = 2 
(figure 7(a)) and d = 3 (figure 7(c)) at sufficiently high temperatures, where 
both the regimes are present. The ratio R of the slopes corresponding to the 
two regimes has got values close to 3 for d = 2 and close to 4 for d = 3, as 
suggested by (|2~4f) . The value of h^sp is obtained from the point of intersection 
of the straight lines fitted to the two regimes. At lower temperatures, however, 
the MD region is absent and the phase diagram here is marked by a single slope 
as shown in figures 7(b) and 7(d). 

Once the pulse is withdrawn, the system relaxes to one of the two equilibrium 
states. The closer one leaves the system to the phase boundary {m w — > 0), larger 
is the relaxation time r^. However, unlike the mean field case, the MC relaxation 
time falls off exponentially with \m w \ away from the phase boundary. Figure 8 
shows the growth of r# as m w — > at a particular T and for a particular At. 
Typical number of MCR used to obtain the data is 400 for L = 40 and 25 for 
L = 400. The best fitted curve through the data points shows the relaxation 
behaviour as follows : 

r K ~ /<T, L)e-^ (T)I "H (25) 

where k(T, L) is a constant depending on temperature and system size and n(T) 
is a constant depending on temperature only. It may be noted from (|25| ) that 
r — > k(T, L) as m w — > 0. Therefore the true divergence at the phase bound- 
ary (where m w = 0) of the relaxation time depends on the nature of k(T,L). 
The inset of figure 8 shows the sharp growth of k(T, L) with the system size. 
The relaxation time tr therefore diverges in the thermodynamic limit (L — ► oo) 
through the constant k. It may be noted that this divergence of tr at the dy- 
namic magnetization-reversal phase boundary occurs even at temperatures far 
below the static critical temperature T c °. 

According to CNT, s c oc \h p \~ d and therefore at any fixed T, stronger fields 
will allow many critical droplets to form and hence the system goes over to the 
MD regime. On the other hand, a weaker field rules out the possibility of more 
than one critical droplet and therefore the system goes over to the SD regime. 
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Figure 9 shows snapshots of the spin configurations at different times in both SD 
and MD regimes. The snapshots at t = t + At corresponds to m w ~ O (10~ 2 ). 
h p > hosp in figure 9(a) and a single large droplet is formed whereas h p < hnsp 
in (b) and many droplets are seen to be formed. It may be noticed from figure 9 
that the boundaries of the droplets are flat with very few kinks on it at t = to+At. 
The probability of growth of a droplet along a flat boundary is very small (only 
25% in case of a square lattice) and hence domain wall movement practically stops 
immediately after the withdrawal of the field. This restricts further nucleation. 
It is then left to very large fluctuations to resume the domain wall movement 
and long time is required for the system to come out of the metastable state and 
subsequently reach the final equilibrium state. Thus the effect of the pulse is to 
initiate the nucleation process and the threshold value of the pulse strength is such 
that within the pulse duration it renders the system with droplets almost without 
any kink in it. This observation justifies the sharp growth of the relaxation time 
at the phase boundary. 

The growth of a length scale at the transition phase boundary can be qual- 
itatively shown from the distribution of domains of reversed spins. We define a 
pseudo-correlation length £ as 



? = E ' R 'f n ' , (26) 



where the radius of gyration R s is defined as B? s = Yh=i \ r i ~ r o| 2 /s, denoting 
the position vector of the ith spin of the domain and tq = Y^t=i r i/ S being the 
centre of mass of the domain. As the transition phase boundary is approached, £ 
is observed to grow with the system size as shown in figure 10. Typical number 
of MCR used for obtaining the data is 10 for L = 1000 and 2000 for L = 50. 
This indicates the divergence of a length scale at the phase boundary in the 
thermodynamic limit. It should be noted, however, that £ is not exactly the 
correlation length of the system | I8| . An estimate for the power law growth of 



the actual correlation length £, as the phase boundary is approached in the MD 
region, will be obtained from the finite size scaling study discussed later in this 
section. 

The order of the magnetization-reversal transition changes with temperature 
and with At even along the same phase boundary. The transition is discontinuous 
all along the low T phase boundary, whereas at higher values of T the nature 
of the transition changes from a continuous to a discontinuous one as one moves 
towards higher values of At. For hp(T) -C hnsp{T), the system is brought to the 
SD regime where the order of the transition is observed to be discontinuous. On 
the other hand continuous transition is observed for hp(T) ^> hnspiT) when the 
system goes over to the MD regime. One can look at the probability distribution 
P(m w ) of m w to determine the order of the phase transition. Figure 11 shows 
the variation of P{m w ) as the phase boundary corresponding to a particular 
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temperature is crossed at two different positions (different At). The data are 
averaged over 500 MCR. The existence of a single peak in (a), which shifts its 
position continuously from +1 to —1 as the phase boundary is crossed, indicates 
the continuous nature of the transition. In (b), however, two peaks of comparable 
strength at positions close to ±mo exist simultaneously. This shows that the 
system can simultaneously reside in both the phases which is a sure indication for 
a discontinuous phase transition. On phase boundaries corresponding to higher 
temperatures the crossover from the discontinuous transition to a continuous one 
is not very sharp and there exists a region around hp = hosp on the phase 
boundary, over which the nature of the transition cannot be determined with 
certainty. This is evident from figure 7, where the data points near the tricritical 
point do not fit to the slope of either of the straight lines corresponding to the 
two different regimes. 

In the region where the transition is continuous in nature one can expect 
scaling arguments to hold. We assume power law behaviour in this regime both 
for m,„ 

P 



m,, 



h, 



n p 



(At, T) 



(27) 



and for the correlation length 



£~ h p -h c p (At,T) 



(28) 



For a finite size system, hp is a function of the system size L. Assuming that at 
the phase boundary £ can at the most reach a value equal to L, one can write 
the finite size scaling form of m w as jl9] : 



m„ 



L-^f [(h p — h^ (At, T, L)) L 1 ^ 



(29) 



where f(x) ~ x^l v as x — > oo. A plot of m w j L~^l u against (h p — h c p (At, T, L)j L 1 ^ 
shows a nice collapse of the data corresponding to L = 50, 100, 200, 400 and 800 
for d = 2 and L = 10, 20, 40, 80 and 120 for d = 3 as shown in figure 12. Typical 
number of MCR used to obtain the data is 5120 for L = 50 in d = 2 and 10000 
for L — 10 in d — 3. The values of the critical exponents obtained from the data 
collapse are (3 = 0.85 ± 0.05 and v = 1.5 ± 0.5 in d = 3 and (3 = 1.00 ± 0.05 and 
v = 2.0 ± 0.5 in d = 2, where h p (At, T) was obtained with an accuracy O (10~ 3 ). 
All attempts to fit similar data to the above finite size scaling form obtained in 
the SD regime failed. 

The accuracy with which hp (At, T) is measured, is very crucial for obtaining 
the critical exponents through finite size scaling. The cumulant method intro- 
duced by Binder et al. is one of the reliable methods which can be employed 

The fourth order cumulant is defined as 



to obtain the value of h c p 



g(L) 



3- 



( m i) " 

( m i) 2 . 



(30) 
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where (m™) = / m™P(m w )dm w . The quantity g(L) is dimensionless and is equal 
to unity for \m w \ ^> 0, while g(L) — > for m w — > 0, assuming a Gaussian distri- 
bution of m m around on the phase boundary Figure 13 shows a plot of g(L) 
against h p at a fixed At and T and the value of the pulse strength correspond- 
ing to the point of intersection of the different curves gives h p (At, T); assuming 



9 = 9 



V 



h„ — hi 



for L = 50 and 2500 



Typical number of MCR used to obtain the data is 50000 

'or L = 800. It is to be noted that none of the curves touch 
the abscissa which corresponds to m w = which is numerically unattainable. 
The closer one gets to m w = better the accuracy in the measurement of h p . In 
principle the minima of g(L) correponding to different L should occur at the same 
position (at h p = h c p ). The shift in the position of the minima of g(L) in figure 
13 is caused by the presence of large fluctuations in measuring higher moments 
of m w . However, this estimate of h p , when used in the scaling fit of (p9|), did not 
significantly improve the estimates of the critical exponents /3 and v. 



4 Summary and Conclusions 

In this paper we have discussed in detail almost all the studies that have been 
made so far on the dynamic magnetization-reversal transition in the Ising model 
under finite duration external magnetic field competing with the existing order 
for T < T c °. Any combination of the pulse strength and duration above the phase 
boundary in the h p — At plane leads to the transition from one ordered phase to 
the equivalent other. We solved numerically the mean field equation of motion 
for the magnetization to obtain the MF phase boundary where the susceptibility 
and the relaxation time were observed to diverge. The divergence of both the 
time {Jr F ) and the length scale (£ MF ) at the MF phase boundary was observed 
even from the analytic solution of the MF equations of motion under a linear 
approximation. Under this approximation, the dynamical critical exponent was 

found to have a value 2 : Tr F ~ ~ — In \m w \, where m w (h p , At, T) = 

gives the phase boundary. The same transition has been studied using Monte 
Carlo simulations in both two and three dimensions. The obtained phase diagram 
is fully consistent with the classical nucleation theory. The nucleation process is 
initiated by the external magnetic field and depending on the strength of the field 
the system nucleates either through the growth of a single droplet or through the 
growth and subsequent coalescence of many droplets. For h p > hasp the system 
belongs to the multi-droplet regime and the transition is continuous in nature; 
whereas for h p < h DSP the system goes over to the single-droplet regime where 
transition is discontinuous. Expecting power law behaviour for both m w and £ in 
multi-droplet regime, the finite size scaling fits give the estimates of the critical 
exponents (3 and v for both d = 2 and 3. Unlike in the MF case, where the 
relaxation time t^ f shows a logarithmic divergence, tr in MC studies falls off 
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exponentially away from m w = and the divergence in tr comes through the 
growth of the prefactor k in ( p5|) with the system size. 

The symmetry breaking transition of the dynamic hysteresis in pure Ising 
systems under oscillating external fields 0@]|f|, where the m — h loop becomes 
asymmetric due to the fact that the magnetization m(t) fails to follow even the 
phase or sign of the rapidly changing field h(t), leads to a dynamic transition. 
This dynamic transition has been studied employing finite size scaling theory 
00 an d the estimates of the critical exponents seem to be consistent with the 



static Ising universality class ||21|| . Although this transition as well as the one dis- 
cussed in this paper occur due to the failure of the system to get out of the 'free 
energy well' corresponding to the existing order because of the lack of proper com- 
bination of the pulse strength and duration, they belong to different universality 
classes. 
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Figure Captions 

Figure 1. Typical time variation of the response magnetizations m(t) for two dif- 
ferent field pulses h(t) with same At and T are shown. The quantities of interest 
to characterize the response magnetizations for both the pulses are indicated. 
Figure 2. MF phase boundaries for three different temperatures. The solid 
line is obtained from numerical solution of (0) and the dotted lines give the 
corresponding analytical estimates in the linear limit. 

Figure 3. Logarithmic divergence of r^ F across the phase boundary for T/T c = 
0.9. The data points shown by circles are obtained from the solution of (0) and 
the solid line corresponds to the solution of the linearized MF equation. 
Figure 4. Divergence of Xg=o across the phase boundary obtained from the 
numerical solution of QT7|). 

Figure 5. Plot of \ q against m w for different values of q. The inset shows the 

linear variation of [i MF ^j against [ln|m w |] _1 . The data points for £ AIF in the 
inset are obtained from the slope of the best fitted straight lines through a plot 
of InXg against q 2 for different values of m w . 

Figure 6. Phase boundaries obtained from the MC study for (a) square lattice 
with L = 100 and (b) simple cubic lattice with L = 50. 

Figure 7. Plot of In At against (h p ) l ~ d along the MC phase boundary, (a) 
T/T c = 0.31 and (b) T/T c = 0.09 for square lattice and (c) T/T c = 0.67 and (d) 
T/T c = 0.11 for simple cubic lattice. The slope ratio R ~ 3.27 in (a) and ~ 3.97 
in (c). 

Figure 8. MC results for the divergence of r R for L = 40, 50, 100, 200 and 400. 
The best fitted straight lines are guide to the eye. The inset shows the variation 
with L of the peak height k in the prefactor of tr in (|25|) . 

Figure 9. Snapshots of spin configurations in a 100 x 100 square lattice at 
different stages (t — t , t x and t + At ) of nucleation, where t < t\ < t + At. 
The dots correspond to +1 spin state, (a) h p = 0.55, At = 300 at T/T c = 0.44 
(SD regime) and (b) h p = 0.52, At = 9 at T/T c = 0.88 (MD regime). 
Figure 10. Variation of £ with L for L = 50, 100, 200, 400, 800 and 1000 for 
MC study on a square lattice. 

Figure 11. Plot of P(m w ) against m w as one crosses the phase boundary for the 
MC study on a 100 x 100 square lattice in (a) MD regime and (b) SD regime. 
Figure 12. Finite size scaling fits : (a) for d — 2 at T/T c = 0.88 and (b) for 
d = 3 at T/T c = 0.67. 

Figure 13. Plot of g(L) against L for L = 50, 100, 200, 400 and 800 in the MD 
regime for MC study on a 100 x 100 square lattice at T = 2.0 for At = 5. 
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